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We perform simulations of Cherenkov Telescope Array (CTA) observations of a young super¬ 
nova remnant RX J1713.7-3946. This target is not only one of the brightest sources ever dis¬ 
covered in very high-energy (VHE) gamma rays but also well observed in other wavebands. In 
X-rays, the emission is dominated by synchrotron radiation, which links directly to the existence 
of high-energy electrons. Radio observations of CO and Hi gas have revealed a highly inhomo¬ 
geneous medium surrounding the SNR, such as clumpy molecular clouds. Therefore gamma rays 
from hadronic interactions are naturally expected. However, the spectrum in GeV energy range 
measured by FermilUPTf indicates more typical of leptonic emission from accelerated electrons. 
Despite lots of multi-wavelength information, the competing interpretations have led to much un¬ 
certainty in the quest of unraveling the true origin of the gamma-ray emission from RX J1713.7- 
3946. CTA will achieve highest performance ever in sensitivity, angular resolution, and energy 
resolution. We estimate CTA capability to examine the emission mechanisms of the gamma rays 
through simulated spatial distribution, spectra, and their time variation. 
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1. Introduction 

More than 100 years have passed since the discovery of cosmic rays but its origin has been 
long in question despite many observational and theoretical researches. The observed spectrum 
of cosmic rays is a power-law shape and has a break around 10*^'^ eV which is so-called “knee”. 
Cosmic rays below the knee energies are thought to be accelerated somewhere in our Galaxy. One 
of the most probable candidates are supernova remnants (SNRs), where the diffusive shock accel¬ 
eration may work at the shock front of SNR blast waves. Evidence for electron acceleration has 
been identified by the detection of synchrotron emission with spatially thin filamentary structures at 
shells of young SNRs (e.g., [1]). On the other hand, gamma rays with hadronic origin was detected 
by Large Area Telescope (EAT) onboard Fermi from middle-aged SNRs 1C443 and W44, which 
are known to be interacting with molecular clouds (MCs)[2]. The observed gamma-ray spectra are 
interpreted as neutral pion decay, which is characterized by a cutoff below 300 MeV, due to the in¬ 
teraction between accelerated cosmic-ray hadrons and MCs. However, it is also observed that that 
the gamma-ray spectra are suppressed above 100 GeV. We hence expect that young SNRs could be 
more plausible as cosmic-ray accelerators to PeV energies, i.e. PeVatron. 

The Cherenkov Telescope Array (CTA) is a next-generation of Imaging Air Cherenkov Tele¬ 
scopes (lACT) observatory which consists of array of the large, middle, and small-sized telescopes 
expanding over km^ area[3, 4]. With higher performance in comparison to the current genera¬ 
tion lACTs, such as better spatial resolution and sensitivity, a search for cosmic-ray PeVatron is 
one of the major scientific objecfives of fhe CTA. A young SNR RX J1713.7—3946 is one of fhe 
brightesf Very High-Energy (VHE) gamma ray sources and spatially exfended emission was ob¬ 
served [5, 6, 7, 8]. The VHE gamma-ray specfrum of RX J1713.7-3946 is fhe mosf precisely 
measured over a wide energy band from 0.3 fo 100 TeV Besides, plenfy of mulfi-wavelengfh ob¬ 
servations have been performed. Term/measured fhe gamma-ray specfrum of RX J1713.7-3946 in 
fhe 3-300 GeV energy range, where fhe observed photon index of 1.5 ± 0.1 is favorable for inverse- 
Compfon emission from accelerafed elecfrons wifh a specfral index of 2.0[9]. X-ray emission is 
dominafed by synchrofron radiafion wifh a good spatial correlafion befween VHE morphology, 
allhough fhe angular resolufion is nof good enough fo be conclusive. On fhe ofher hand, radio 
observations of CO and Hi gas have revealed a clumpy molecular clouds (MCs) surrounding fhe 
SNR, and reported evidence for interaction befween fhe MCs and fhe SNR shock (e.g. [10, 11, 12]. 
The hadronic gamma-ray emission is nalurally expecled to reproduce fhe obseved spelrum (e.g. 
[13]) . However, our idea here is lhal, if fhe hadronic gamma rays do exisf, such componenl mighl 
be hidden by fhe dominanf leplonic gamma-ray emission. If is of a greaf inleresl whelher fhe im¬ 
proved sensitivity of fhe CTA could deled fhe possible bul dim hadronic gamma rays. Hence Ibis 
objecl is a very good fargel for deep observations, and also for conslraining Iheorelical models of 
cosmic-ray acceleralion. 

2. Aims and methods of simulations 

The major purpose of our simulalion studies is to show an example of analysis strategy when 
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we will obtain real data, and to evaluate the capabilities of CTA on finding a clue for the hadronic 
gamma rays. 

First we perform morphological analyses in order to find ouf fhe dominanf componenf of fhe 
VHE gamma-ray emission from RX J1713.7-3946. As for fhe hadronic gamma rays, fhe morphol¬ 
ogy should be relafed fo spatial disfribufion of fhe accelerafed protons and fhaf of fhe inferacfing 
mailer densify which is indicaled by fhe CO and Hi morphology obfained by fhe radio observations. 
Since we currenfly do nof know fhe CR disfribufion, we here roughly assume lhal CR would be 
filled homogeneously inside fhe SNR. On fhe olher hand, fhe morphology of fhe leplonic gamma- 
ray emission may be fraced by fhaf of synchrofron X-ray. We should nofe fhaf fhe X-ray and VHE 
morphologies are nof always complefely same. The brighfness of fhe synchrofron emission reflecfs 
nof only for spalial election dislribution bul also local magnetic fields. And fhe inverse-Compfon 
process is also coupled fo energy densities of large! phofons which includes infrared and/or oplical 
pholon field in addition fo fhe cosmic microwave background. However, the overall structure could 
be approximated by the X-ray morphology. We hence apply the radio and X-ray images as tem¬ 
plates for the hadronic and leptonic gamma-ray morphology, respectively. We perform maximum 
likelihood test to quantitatively determine which component dominates the VHE emission from the 
SNR. The templates also contain spectral information. Here we simply assume the same spectral 
shape is assumed over the SNR image. The spectrum of the leptonic component is modeled as 



(2.1) 


where Ai is a normalization factor. Eg is a photon index, and is a cutoff energy. Input values 
for Eg and are 2.04 and 17.9 TeV, respectively, as reported by H.E.S.S. observations[8]. For the 
hadronic emission, the spectrum is described as follows. 



( 2 . 2 ) 


where A 2 is a normalization factor, Tp is a photon index, and Ec is a cutoff energy. We adopt 
Yp = 2.0 and E^ = 300 TeV as fiducial parameters. Therefore Ai and A 2 , or their ratio A 2 /A 1 , are 
the parameter to be investigated, requiring that the sum of the integral fluxes befween 1 and 10 TeV 
are equal fo thaf measured by H.E.S.S.. If fhe hadronic gamma ray is greafer, we could conclude 
fhaf dominanf par! of fhe hadronic componenf were nof accelerafed fo fhe knee energies in RX 
J1713.7-3946. Searching for a spectial componenf lhal extends to PeV energies, we also look for 
a dimmer hadronic component by spectral analysis. The maximum likelihood fit will be performed 
in order to unfold spectra for each component and evaluate statistical significance of fhe hadronic 
gamma-ray defeclion. Here fhe spatial templates are also considered to calculale fhe likelihood. 

We also evaluate fhe capabilily of delecling fhe lime varialion of fhe spectial culoff energy, 
Emax! wilh longer time scale. The maximum energy of fhe CR speclrum is delermined by a balance 
among acceleration, cooling and escape. Hence Emax variation, increase or decrease, depends on 
fhe SNR age and also on acceleration Iheories. In fhe case of RX J1713.7-3946, ~10% variation 
in 10-20 years may be expecled, where E^ax could vary faster in fhe leptonic scenario[14]. Since 
fhe CTA will be operational for a few tens of years, such a long-lerm sludy in VHE energies will 
become possible. This may be a unique approach for idenlifying fhe VHE emission mechanism. 
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(a) CTA leptonic dominant case (A2/Al=0.01) 
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(b) CTA hadronic dominant case (A2/Al=100) 
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Figure 1: Simulated gamma-ray images of (a) A 2 /A 1 = 0.01 (leptonic dominant case) and (b) A 2 /A 1 = 100 
(hadronic dominant case) with Fp = 2.0 and El’ = 300 TeV. The green contours show (a) XMM-Newton 
X-ray intensity [16] and (b) total interstellar proton column density [11], which smoothed to match the PSF 
of CTA. The unit of color axis is counts pixel^* for both panels. 


The simulation software package that we use in this study is ctools version 00-07-01 [15]. We 
use a preliminary instrumental response function which corresponds to the CTA southern array lo¬ 
cated at the candidate site Aar (southern Namibia). When we perform the simulations, the Galactic 
diffuse emission and isotropic background due to gamma-like charged cosmic rays are taken into 
account for background photons in the field of view. 


3. Results 
3.1 7 -ray image 

In order to clarify the imaging capability of CTA, we first intend to simulate different gamma 
ray images in the energy range of 1-100 TeV by tuning A 2 /A 1 , the ratio between the hadronic 
and leptonic gamma rays. Figures la and lb show the simulated gamma-ray images in leptonic 
dominant case (A 2 /A 1 = 0.01) and hadronic dominant case (A 2 /A 1 = 100), respectively. Each 
gamma-ray image is similar to each overlaid contour, which corresponds to the non-thermal X 
rays and the total ISM protons including both molecular and atomic hydrogen, respectively. On 
the other hand, the spatial distributions of gamma rays are apparently different from each other, 
particularly the north and the southwest. In this extreme case, we can therefore determine the 
major component of the VHE emission by the morphological study with CTA. Incidentally, we 
found that A 2 /A 1 = 1-10 showed the best spatial correspondence with the H.E.S.S. excess counts 
map [8] with a correlation coefficient of ~0.7-0.8. We continue to study the systematic error 
estimation and quantitative evaluation for the morphological difference. 
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Figure 2: Spectral energy distribution of the gamma ray emission obtained by analyzing the CTA simulation 
data for RX J1713.7—3946 with A 2 /A 1 = 0.1. The blue and red squares are the spectral points for the 
leptonic and hadronic spatial templates, respectively. Only statistical errors are presented. The black squares 
are the total fluxes of the leptonic and hadronic components. The black vertical bars are the errors for the 
total fluxes obtained by adding the errors for two components in quadrature. The blue, red, and black solid 
lines show the input spectra for the leptonic component, the hadronic component, and the total, respectively. 


3.2 Spectrum 

Assuming that the leptonic component is dominant, we subsequently proceed to search for 
a “hidden” hard component with a hadronic origin. Using reconstructed energy band of 0.5 - 
100 TeV with 50-hour observation, likelihood analyses show the significance > 10a to observe 
a dimmer hard component even for a small A 2 /A 1 = 0.02. Note that the result may be rather 
optimistic since the fitting templates are the same as the input for the simulation. 

We then proceed to perform maximum likelihood fittings for the simulation data (with a ratio 
^ 2 Mi = 0.1, for the safety) in 12 logarithmically spaced energy bands. Figure 2 shows the resulting 
spectrum from our ‘bin-by-bin’ analysis of the same 50 hr of simulation data. It is clear that our 
likelihood fits reproduce the simulated spectrum for each spatial template (i.e. hadronic or leptonic 
morphology), which demonstrates the capability of detecting the hidden hadronic component in 
the best case scenario. 

3.3 Time variation of cutoff energy 

Detecting the time variation of the cutoff energy of the gamma ray spectrum can provide a 
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Figure 3: Significance of the detected variation of Eq as a function of exposure time. Dashed lines represent 
the best-ht curve to each dataset which is proportional to v^. Squares and circles represent results for the 
AEc/Ec = —10% and +10% case, respectively. 


clue to the emission scenario and acceleration theories. The simulations are performed for three 
sets of intrinsic cutoff energies at 17.9 (nominal), 19.7 (+10% case) and 16.1 TeV (—10% case) 
each. As a start, we consider a variation of AEc/Ec = ±10% to show how sensitive CTA will be 
to such fractional changes in the spectral cutoff. Here we consider the pure leptonic scenario as an 
example, and used 0.2-100 TeV photons. 

We define a significance, s, for the observed En variation as 5±(t) = where En is 

the nominal value and E± is the best fit Ec for the cases with a ±10% variation. Gq and a± are the 
corresponding errors. We repeat all of our simulations for 100 times and take the average of the 
calculated s{t) for each run, and show the result on Figure 3 . Our result indicates that a decrease of 
Ec is slightly easier to identify than an increase. As a result, a lower cutoff energy can actually be 
easier to measure and precisely for a given exposure. If we observe > 60 hrs in the two epochs, we 
are able to achieve a 3a detection for the AEc/Ec = —10% case, whereas ~ 70 hrs are necessary 
for the +10% case. 

4. Summary 

In this paper, we have briefly introduced our feasibility studies for CTA observations of RX J1713.7- 
3946, mostly with 50 hrs observations. We showed that a 50-hr observation may be enough to 
identify the dominant gamma ray emission component by the morphology obtained with CTA. 

And in the case that the leptonic emission would be dominant, we should be able to quantify both 
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the leptonic and hadronic components through spectral analysis if they are mixed with a ratio of 
^ 2/^1 =0.1 or less. Interestingly, we also found that CTA will be able to reveal variations of the 
spectral cutoff energy over 10-20 years, for the very first time. A variation of AEc/Ec = ±10% 
could be detected provided that an exposure time longer than 70 hr' can be secured for the two 
epoch. 

However we know our present study is based on a fairly simplified input model and may 
contain systematic errors of which estimation is not trivial. And also this study can be extended to 
use more theoretically justified models. More studies are left for future works. 
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